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ABSTRACT 



The influence of the uncertainties in the rate coefficient data for electron-impact excitation and ionization on non-LTE 
Li line formation in cool stellar atmospheres is investigated. We examine the electron collision data used in previous 
non-LTE calculations and compare them to recent calculations that use convergent close-coupling (CCC) techniques 
and to our own calculations using the R-matrix with pseudostates (RMPS) method. We find excellent agreement 
between rate coefficients from the CCC and RMPS calculations, and reasonable agreement between these data and 
the semi-empirical data used in non-LTE calculations up to now. The results of non-LTE calculations using the old 
and new data sets are compared and only small differences found: about 0.01 dex (~ 2%) or less in the abundance 
corrections. We therefore conclude that the influence on non-LTE calculations of uncertainties in the electron collision 
data is negligible. Indeed, together with the collision data for the charge exchange process Li(3s) + H ^ Li + + H~ now 
available, and barring the existence of an unknown important collisional process, the collisional data in general is not 
a source of significant uncertainty in non-LTE Li line formation calculations. 
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1. Introduction 

Lithium abundances in stellar atmospheres are key ob- 
servables in astrophysics, giving crucial information on 
stellar evolution and mixing, stellar and Big Bang nu- 
cleosynthesis, cosmic ray spallation and perhaps even - 
planet formation (e.g . lLambertlll993t ICarlsson et al.lll994t 
iMelendez et al.l l2010h . Such abundances are interpreted 
from observations of the few Li I lines found in stellar spec- 
tra, and to obtain accurate results, it is important to ac- 
count f or departures from local thermodynamic equilibrium 
(LTEi (ISteenbock fc Holwegerl [l98l ICarlsson et aTlfl99i 
iLind et al.ll2009f h Reliable modelling of line formation in 
non-LTE requires detailed knowledge of all important ra- 
diative and collisional processes on the atom of interest, 
Li. 

The collisional processes are particularly challeng- 
in g. For many years fo l lowin g the pioneering study 
of ISteenbock fc Holwegerl (|1984D . the importance of in- 
clastic hydrogen collisions was a major uncertainty. 
Detail ed quantum scattering ca lculations of Li+H colli- 
sions (|Belvaev fc Barklem! 12003]) and application to non- 
LTE modelling (jBarklem et alJl2003t iLiiid et al.ll2009h have 
allowed this question to be answered, and it was found that 
direct excitations by hydrogen collisions essentially have 
no influence. However, a related charge-exchange process 
Li(3s) + H ^ Li + + H~ was shown to be important, result- 
ing in differences in derived abundances of about 0.05 dex 
in solar- metallicity stars and 0.1 dex in metal-poor stars. 
Moreover, it was found that the results were not sensitive to 
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uncertainties in the data for this process. When the rate co- 
efficient was altered by factors in keeping with the expected 
uncertainty in the theoretical calculation, the effects on the 
line formation were practically negligible (< 0.01 dex in 
derived abundances) . 

Thus, barring there being additional important colli- 
sional processes on Li we are unaware of, the uncertainties 
regarding collisional data lie, perhaps, with the electron 
collisions. The question of quantifyi ng these has on l y been 
touched upon briefly in the past. ICarlsson et al.l (|1994[) 
made calculations exploring the sensitivity of their results 
to various input parameters, including collision cross sec- 
tions. They identified the oscillator strengths of the lines of 
interest (6708 and 6104 A) and the photoionization cross 
sections from 2s and 2p as the most important atomic data 
in determining the uncertainties of their calculations. They 
estimate an error in the abundance corrections of less than 
0.01 dex arising from these sources. That the sensitivity 
to collisional cross sections was explored, though not dis- 
cussed explicitly implies even smaller uncertainties due to 
collision cross sections, but de tails have not been published . 
This study and others (e.g. ISteenbock fc Holwegerl Il984t 
ILind et aLl I2009T) make use of em piric ally correcte d calcu- 
lations by Ivan Regemorteil ()1962l ) and iParkl ()1971l ) for ex- 
citation and lAllenl (|1976lf or similar for ionization. 

Given the great astrophysical importance of Li abun- 
dances and the considerable advances in calculation meth- 
ods for electron scattering in the intervening period, in par- 
ticular advanced close-coupling methods that are able to 
account for effects of coupling to the target continuum, we 
considered it worthwhile and prudent to accurately quan- 
tify the uncertainties associated with the electron collision 
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data. In this paper we examine the electron collision data 
in the literature, and perform a new i?-matrix calculation, 
in order to accurately estimate the uncertainties in such 
calculations. The data are then used in non-LTE Li I line 
formation modelling to assess the resulting uncertainties in 
stellar Li abundances. 

2. Electron collision data 

In this section we examine the existing data for excitation 
and ionization by electron impact used in non-LTE cal- 
culations. We compare with advanced close-coupling cal- 
culations including our own calculations for the excitation 
processes and present the details. 

2.1. Excitation 

There are a number of existing calculations for electron- 
impact excitation Li(ni) + e — »■ Li(n'Z') + e, though only 
relatively recently have there been any large-scale calcula- 
tions including a significant number of excited states. Two 
advanced close-coupling methods have been used, namely 
convergent close-coupling (CCC, iBrav fc Stelbovicsl 19921) 
and i?-matrix with pseudostates (RMPS. iBartschat et al.l 
1996) approaches. The advantage of the CCC method is 
that, at least for collision systems with two electrons, it 
has been shown to yield convergent results as basis size is 
increased; however, it has the disadvantage of requiring con- 
siderable computational effort. The advantage of the RMPS 
method is that the i?-matrix is calculated once and then re- 
sults for many collision energies can be calculated with lit- 
tle additional effort. This is particularly important for the 
study of resonances requiring dense energy grids. The dis- 
advantage, however, is that convergence with basis size can 
be slower resulting in the ap pearance of pseudo- resonances. 

iSchweinzer et al.1 (|1999D have published data from ex- 
tensive 45-state CCC calculations for transitions with n < 3 
and n' < 4, and include semi-empirical cross sections for 
n = 4, n' = 4. The calculations show good agreement with 
experiments where availabl e, and analyt i c fits to the cross 
section data are provided. iGriffin et al.l (|200ll ) have done 
calculations using a 55-state RMPS approach resulting in 
cross section data for all transitions with n — 2, I = and 
n' < 4. Their calculations demonstrate the importance of 
including the coupling to the target continuum via pseu- 
dostates, particularly at collisions energies greater than the 
ionization energy, i.e. > 5.4 eV, where the effects are small 
for 2s —> 2p, but become significant for 2s — > 3/ and even 
greater for 2s — > Al. The cross sections are in good agree- 
ment with the CCC and experimental results where com- 
pared. Since there is a possibility that resonances at low 
collision energies are unresolved by the CCC calculations 
and could contribute significantly to the rate coefficients, 
we performed our own RMPS calculations for n < 4 and 
n' < 4 and n'V — 5s, which gives the additional advantage 
of an independent check. These calculations are described 
below in Sect. 12.1.11 and compared with the existing data 
in Sect. l2"7L~2l 

A s men tioned above, the rate coefficients calculated by 

iParkl ([1971D have been used in most calculations of non-LTE 
formation of Li I lines in cool stars (|Steenbock fc Holwegei] 
1984tlCarlsson et al.lll994l:lLind et al.ll2009USbordone et all 
20iq iTakeda fc Kawanomotol 12005^ 7" so they repre- 
sent an important basis for comparison in this work. 



Table 1. Comparison with experiment for the energies of 
the lowest states of Li I. 

State Energy (Ry) 

This Work Experiment Difference 



Is 2 2s 


2 S 








2p 


2po 


0.1350 


0.1358 


-0.59% 


3s 


2 S 


0.2474 


0.2479 


-0.20% 


3p 


2po 


0.2807 


0.2818 


-0.39% 


3d 


2 D 


0.2839 


0.2851 


-0.42% 


4s 


2 S 


0.3180 


0.3191 


-0.34% 


4p 


2po 


0.3312 


0.3323 


-0.33% 


id 


2 D 


0.3325 


0.3337 


-0.36% 


4/ 


2po 


0.3326 


0.3338 


-0.36% 


5s 


2 S 


0.3479 


0.3490 


-0.32% 



Notes. Energies are relative to the ground state and compared 
with experimental values obtained from NIST data (see text). 
The energy difference is measured (-E-rhcor — £Expcr)/£Wpcr. 



ISteenbock k, Holweeerl (|l984l) replaced th e data of the im- 
portan t 2 s — > 2p transition with that of Ivan Reeemorterl 
(|1962l ). Both these sources are based on Born approxima- 
tion calculations, which are well known to be valid only at 
high impact energies and to substantially overestimate the 
cross sections near threshold. Thus, empirical corrections 
determined from various experiments are applied. Based 
on comparison with the experiments, both works conclude 
that the calculated rates are not likely to be in error by 
more than a factor of two. 

2.1.1. R-matrix calculation 

We made a 34-state RMPS calculation of the excitation by 
electron impacts, using existing freely available computer 
codes. These calculations will now be described. 

For the ini tial atomic structur e calculations we used the 
cod e CIVPOL dLe DourneuJ Il976l : iHibbert et"aT1 Il977t see 
also [P lummcr ~etalll2004l) . which is an adaptation of CIV3 
that allows the construction of polarization 
pseudostates. Such pseudostates are used to account for 
coupling to the target continuum and are needed to ob- 
tain a correct description of the dipole polarization and 
thus the l ong-range interactio n potential. They have been 
shown by IGriffin et all (|2001l ) to be important at interme- 
diate energies. CIV3 and CIVPOL use the configuration in- 
teraction (CI) method and the radial parts of the (pseudo) 
orb itals are repres ented by Slater- type orbita ls; for details 
see IHibbert! (| 19751) and IHibbert et all (|l977l) . The atomic 
structure of the target Li I was built by optimising spectro- 
scopic orbitals with all allowed nl values up to nl = 5/; 
the Is and 2s orbitals were H artree-Fock orbitals taken 
from IClementi fc Roettil (|l974l ). The energies obtained for 
the ten lowest-lying spectroscopic states are compared with 
experimen tal values ta ken from the NIST Atomic Spectra 
Database (jNISTl 120101 ) in Table [TJ and the agreement is 
satisfactory. 

Once the spectroscopic orbitals (and states) were calcu- 
lated, we added n — 6(s,p,d) and 7{s,p,d) pseudo-orbitals 
to describe the polarization of the ground and first excited 
state, respectively. To test the quality of the polarization 
pseudostates, we compared the static dipole polarizabilities 
obtained with experimental and theoretical values found 
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Table 2. Dipolc polarizabilitics of the ground and first ex- 
cited state compared with values found in the literature. 



State 


Polarizability (ao) 




This Work Other Works 


Is 2 2s fS) 


168.5 (a) 162.3, (b) 164.0(34) 




(c) 164.1 


Is 2 2p ( Z P C ) 


148.5 (a) 117.8 



Ashbv & vai^W^ngaardenl (|2003f ) using 
' Molof et all (11974) usin g E-H gra- 



Pipin fc Bishod (jl992l ) using CI 



References, (a 

Coulomb Approximation; (b) 
dient balance technique; (c) 
Hylleraas. Bold values are experimental. 



in the literature, see Table [U For the polarization of the 
ground state there is good agreement between our results 
and those from the literature. In the case of the polariz- 
ability of the first excited state, we found only on e other 
theoretical result (| Ashbv fc van Wiingaa rdcn 2003|, where 
the polarizability is calculated in the Coulomb approxima- 
tion. 

The computed orbitals were then used in a RMPS calcu- 
lation of the ele ctron-impact excitation cross sections. The 
RMATRX I code (jBerrineton et al.lll995h was used for the in- 
ternal region problem, and for the extern al region problem 
we em ployed a version of the code STGF (|Berrington et all 
1 19871) modified for treating collisions with neutral atoms 
(|BadnellllT999l) . 

After exploration of the dependence of our results on 
the calculation details, we finally adopted the following. 
The target atom consists of 30 spectroscopic states (the 
ten states in Table [1] with an additional 20 spectroscopic 
states) and four polarization pseudostates; 34 states in to- 
tal. Target atoms with 40 and 45 states were tested and 
showed no significant change in the cross section results. 
For both the iV-elcctron and N + 1-electron configurations, 
maximum occupation numbers of 2 for the spectroscopic or- 
bitals and 1 for the pseudo-orbitals were adopted, with dou- 
ble excitations permitted. Calculations extending the N+l- 
electron configurations to allow all possible ways of adding 
one electron to the A^-electron configurations showed no 
significant differences in the inelastic cross sections. The R- 
matrix boundary between the inner and outer regions was 
set at 83.4 a.u. Partial waves up to L = 50 were calculated. 

2.1.2. Results and comparison 

In Fig.[T]some of the resulting cross sections are shown, and 
then compared with a calculation with out pseudostates, 
and w ith the fits to the CCC results of ISchweinzer et "al] 
(|l999h . It is seen that the final cross sections compare well 
with the CCC results, particularly at low collision energies, 
of greatest interest for the rates at cool star temperatures. 
It appears that our results are not always fully converged 
at intermediate energies. For example, at collision energies 
above 5.4 eV, where processes involving the continuum are 
directly accessible, the 2s — > 3p cross section lies between 
the CCC result and the calculation without pseudostates. 
Griff in et al.l (|2001l) obtain good agreement between their 
somewhat larger 55-state RMPS calculation and the CCC 
result. Thus, we suspect that the intermediate-energy re- 
sults are not converged in this and in some other transitions; 
we note that this may be a result of the smaller basis size, 



but could also come from other details of the calculation. 
However, as we will see, this is not important for calculat- 
ing the rate coefficients at the temperatures and precision of 
interest for stellar atmosphere applications. The behaviour 
within a couple of eV of threshold is far more important. 
Examining Fig.[T]we find generally good agreement between 
the CCC results and ours in this region. 

Rate coefficients were computed both from our 
RMPS cross sec t ions a nd the CCC cross sections from 
ISchweinzer et al.l (| 19991) by integrating the cross sections 
folded with a Maxwellian velocity distribution. The rate 
coefficients from our RMPS calculations are presented in 
Table [3] for a range of temperatures. In the table the ex- 
citation rate coefficients are presented and de-excitation 
rates may be obtained from the detailed balance relation. 
The two sets of rates are compared in the leftmost panel 
of Fig. [21 and the agreement is excellent. This suggests 
that neither the inclusion of resonances nor the proper 
convergence at intermediate energies is vital at tempera- 
tures of interest, though we note that some cancellation 
may occur. Assuming a log-normal distribution for the ratio 
(cw)this work/(cw)ccc we find a mean of 0.83 and a scatter 
of 1.86; i.e. an offset of less than 20% and a scatter of less 
than a factor of 2. There is a tendency for our rate coef- 
ficients to be lower than the CCC ones for the transitions 
with the highest rate coefficients. Of these seven transitions 
grouped with rate coefficients greater than 10 -6 cm 3 s _1 , 
the one lying roughly on the one-to-one relation is 3p — > 3d. 
The remaining six transitions correspond to the n = 4 
and n' = 4 transitions , where, as mentioned earlier, the 
ISchweinzer et all (|l999l ) data are not actually CCC calcula- 
tions, but semi-empirical data. The agreement for the most 
important low-lying transitions (shown as different sym- 
bols) is excellent. 

Ou r data are also compared with the data from iParkl 
(| 19711 ) and data calculated u sing the van Regemorter for- 
mula (|van Regemorterl Il962l ). the comparisons shown in 
the middle and rightmost panels of Fig. [51 respectively. 
Comparison with the data of Park is generally good, al- 
though the scatter is significantly larger than seen for the 
CCC rates. The Park values are typically around 40% 
greater than ours, and the scatter is around a factor of 
6. It is also notable that the rate for 2p — s> 3s, correspond- 
ing to the 8126 A line, is more than 3 orders of magni- 
tude greater in our calculations than found by Park. The 
4c? — > 4/ is the other significantly discrepant transition. The 
van Regemorter formula is only valid for optically allowed 
transitions, showing a smaller number of comparisons. The 
comparison with van Regemorter has an even greater scat- 
ter than found for Park, around a factor of 18, though owing 
to the small number of comparisons the scatter is domi- 
nated by a small number of outliers. The offset is slightly 
smaller than for the Park data. Van Regemorter's formula 
gives a rate for the 2p — > 3s transition in much better agree- 
ment with the rate from our calculation and the CCC data, 
while a few other transitions show large discrepancies. 

The non-LTE line source function for the very impor- 
tant 6708 A res onance line is describ ed well by a two-level 
approximation (jCarlsson et al.| [l994) , and thus the data for 
2s — » 2p are of special interest. The data are compared in 
Table [4] and all f our ca l culati ons are in reasonable agree- 
ment, though the IParkl (|1971l) value is roughly a factor of 
two larger than the other three. 
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Table 3. Rate coefficients for inelastic collisions between electrons and Li I atoms for different temperatures. 



Initital 














Final state 














state 


2p ( 2 P°) 


3s ( 2 S°) 


3p ( 2 P°) 


3d ( 2 D°) 


4s ( 2 S°) 


4p ( 2 P°) 


4d ( 2 D e ) 


4f ( 2 F°) 


5s ( 2 S°) 
















1 000 K 














2s ( 2 S°) 


2.04(-16) 


1.93(- 


-24) 


7.55(- 


27) 


3.51(-27) 


1.74(-29) 


6.23(-31) 


4.79(-31) 


4.76( 


-3D 


3.18 


— 32) 


2p ( 2 P°) 




8.91(- 


-15) 


5.30(- 


17) 


4.72(-17) 


2.81(-20) 


4.34(-21) 


3.17(-21) 


3.08 


-21) 


7.00 


-23) 


3s ( 2 S°) 








6.34(- 


09) 


3.94(-09) 


7.58(-12) 


7.02(-13) 


5.47(-13) 


7.16 


-13) 


2.32 


-14) 


3p ( 2 P°) 












1.12(-05) 


1.42(-09) 


2.07(-10) 


1.47(-10) 


1.63 


-10) 


2.05 


-12) 


3d ( 2 D°) 












_ _ 


1.62(-09) 


2.59(-10) 


1.96(-10) 


2.76 


-10) 


3.10 


-12) 


4s ( 2 S°) 














_ _ 


2.14(-07) 


1.88(-07) 


2.91 


-07) 


5.41 


-09) 


4p ( 2 P°) 
















_ _ 


3.79(-05) 


1.33 


-05) 


3.58 


-08) 


4d ( 2 D°) 





















6.11 


-05) 


4.47 


-08) 


4f ( 2 F°) 














2 000 K 










4 57 


— 08) 


2s ( 2 S°) 


8.54(-12) 


4.78(- 


-16) 


2 17(- 


17) 


1.57(-17) 


7.14(-19) 


9.99(-20) 


9.18(-20) 


8.61 


-20) 


1.95 


-20) 


2p ( 2 P°) 




4.17(- 


-11) 


3.13(- 


12) 


3.61(-12) 


2.70(-14) 


1.49(-14) 


1.32(-14) 


1.22 


-14) 


9.37 


-16) 


3s ( 2 S°) 








8.29(- 


08) 


7.68(-08) 


1.40(-09) 


3.92(-10) 


3.84(-10) 


5.52 


-10) 


5.12 


-11) 


3p ( 2 P°) 












1.67(-05) 


1.96(-08) 


9.20(-09) 


7.71(-09) 


9.31 


-09) 


3.56 


-10) 


3d ( 2 D°) 















1.96(-08) 


8.70(-09) 


9.01(-09) 


1.29 


-08) 


4.56 


-10) 

1 


4s ( 2 S°) 














_ _ 


7.34(-07) 


8.27(-07) 


9.39( 


-07) 


4.94 


-08) 


4p ( 2 P°) 
















_ _ 


5.13(-05) 


1.40 


-05) 


1.54 


-07) 


4d ( 2 D°) 





















5.99 


-05) 


1.61 


-07) 


4f ( 2 F°) 














5 000 K 










1 40 


Ul ) 


2s ( 2 S°) 


5.91(-09) 


3.88(- 


-11) 


8.32(- 


12) 


8.06(-12) 


1.16(-12) 


4.52(-13) 


5.23(-13) 


4.24( 


-13) 


2.25 


-13) 


2p ( 2 P°) 




4.97(- 


-09) 


1.91(- 


09) 


2.62(-09) 


7.82(-ll) 


1.00(-10) 


1.18(-10) 


1.06 


-10) 


1.88 


-11) 


3s ( 2 S°) 








4.96(- 


07) 


4.75(-07) 


2.32(-08) 


1.50(-08) 


1.53(-08) 


3.14 


-08) 


4.88 


-09) 


3p ( 2 P°) 












1.96(-05) 


9.03(-08) 


8.17(-08) 


8.09(-08) 


1.24 


-07) 


8.32( 


-09) 


3d ( 2 D e ) 












_ _ 


7.26(-08) 


5.63(-08) 


9.24(-08) 


1.46 


-07) 


8.44 


-09) 


4s ( 2 S°) 














_ _ 


2.71(-06) 


2.47(-06) 


1.65 


-06) 


1.74 


-07) 


4p ( 2 P°) 
















_ _ 


5.32(-05) 


1.14 


-05) 


5.00( 


-07) 


4d ( 2 D°) 





















4.70 


-05) 


3.54 


-07) 


4f ( 2 F°) 














8 000 K 










2 20 




2s ( 2 S°) 


3.25(-08) 


5.44(- 


-10) 


1.86(- 


10) 


2.07(-10) 


3.71(-11) 


2.05(-ll) 


2.80(-ll) 


2.02( 


-11) 


1.25 


-11) 


2p ( 2 P°) 




1.43(- 


-08) 


8.63(- 


09) 


1.32(-08) 


5.39(-10) 


8.85(-10) 


1.24(-09) 


1.11 


-09) 


2.15 


-10) 

1 


3s ( 2 S°) 








9.63(- 


07) 


7.18(-07) 


4.21(-08) 


3.76(-08) 


3.70(-08) 


9.52 


-08) 


1.46 


-08) 


3p ( 2 P°) 












1.86(-05) 


1.43(-07) 


1.42(-07) 


1.57(-07) 


2.75 


-07) 


1.84 


-08) 


3d ( 2 D e ) 















9.27(-08) 


8.23(-08) 


1.72(-07) 


3.30 


-07) 


1.57 


-08) 


4s ( 2 S°) 














_ _ 


4.80(-06) 


3.37(-06) 


1.73 


-06) 


2.43 


-07) 


4p ( 2 P°) 
















_ _ 


4.80(-05) 


9.41 


-06) 


7.60 


-07) 


4d ( 2 D°) 





















3.90 


-05) 


4.05 


-07) 


4f ( 2 F°) 














12 000 K 










2.13 


— 07i 
U i j 


2s ( 2 S°) 


8.67(-08) 


2.15(- 


-09) 


9.95(- 


10) 


1.29(-09) 


2.54(-10) 


1.74(-10) 


2.75(-10) 


1.73 


-10) 


1.12 


-10) 


2p ( 2 P°) 




2.46(- 


-08) 


1.92(- 


08) 


3.36(-08) 


1.58(-09) 


3.06(-09) 


4.98(-09) 


4.19 


-09) 


7.83 


-10) 


3s ( 2 S C ) 








1.56(- 


06) 


8.72(-07) 


5.99(-08) 


6.46(-08) 


6.32(-08) 


1.81 


-07) 


2.62 


-08) 


3p ( 2 P°) 












1.68(-05) 


2.01(-07) 


2.04(-07) 


2.44(-07) 


4.56 


-07) 


2.86 


-08) 


3d ( 2 D°) 












_ _ 


1.01(-07) 


9.88(-08) 


2.49(-07) 


5.94 


-07) 


2.02 


-08) 


4s ( 2 S°) 

















7.10(-06) 


4.04(-06) 


1.69 


-06) 


3.08 


-07) 


4p ( 2 P°) 
















_ _ 


4.19(-05) 


7.77 


-06) 


9.94 


-07) 


4d ( 2 D°) 





















3.24 


-05) 


4.12 


-07) 


4f ( 2 F°) 














20 000 K 










1 90 
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Fig. 1. Comparison of electron- impact excitation cross sections with collision energy for some transitions involving the 
lowest states of Li I. Solid lines show our calculation using the polarization pseudo-orbitals, dashed l ines are the cross 
sectio ns obtained without using pseudo-orbitals, and dotted lines are the fits to the CCC calculations (jSchweinzer et al.l 
1999). 




10" icr' u 10" 1Q ,10 , in" 

Rate Coefficients [cm /sj (CCC) 




10"' 2 10 



1 p? 10" ,1 0" , 1 0" 
Rate Coefficients [cm /sj (Park) 




10" " 10"'° 10"" JO" 6 10" 4 10" 
Rate Coefficients [cm /sj (Van Regemorter) 



Fig. 2. Rate coefficien ts at 8000 K o btained by us compared with those calculated from the CCC calculations 
Sch weinzer et al.l 19991 left), those of (|Parkl [l97ll middle) and those calculated using the van Regemorter formula 



van Regemorterlll962f ). The latter is only valid for optically allowed transitions, hence the limited number of compar- 



ison points. The different symbols are for the rate coefficients corresponding to the transitions: 2s —t 2p [ O, 6708 A ], 
2p -> 3d [ □, 6104 A ], and 2p ->• 3s [ A, 8126 A ]. 



2.2. Ionization 

Regarding electron-imp act ionization Li(n/) + e — > Li + + 2e, 
ISchweinzer et al.l (|l999h have also calculated cross sections 
via the CCC method for the 2p and n = 3 states. For 
ionization cr oss sections from the ground state (2s) cal- 
culations by iLennon et all |l988) are recommended by 



ISchweinzer et~afl (Il999h . The results of CCC and RMPS 
calculations for this proce ss where n = 2 hav e been shown 
to be in good agreement (jColgan et al.ll200ll) . Thus, since 
this process is likely to be less important in this application, 
we adopted the CCC data directly. As mentioned, non-LTE 
calculations have typically used a semi-empirical formula 
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Table 4. Comparison of the rate coefficient (av) for the 
transition 2s -> 2p at 8000 K. 



Method 



(av) (cm 3 



This Work 3.28(-08) 

Schweinzer et al. ( 1999) 3.33(-08) 

Park (1971) 6.36(-08) 

van RcEcmortcr (1962) 2.46(-08) 



from Allen J 1976ft based on a fit to experimental data (see 
lSeatonlll9620 . A comparison of c ross sections from t he sem i- 
empirical formula with those of ISchweinzer et all (|1999ft is 
presented in Fig. [3] The comparison shows a tendency for 
the ionization cross sections of the semi-empirical formula 
to be significantly larger than the CCC results. Rate coeffi- 
cients calculated by integration over a Maxwellian velocity 
distribution are presented in Table [5] 

3. Application to non-LTE Li I line formation 

The most commonly observed lines for Li I in stellar spectra 
are those corresponding to the 2s- 2p and the 2p-3d transi- 
tions at 6708 A and 6104 A, re spectively. In Li-r i ch stars the 
8126 A line is also observed (jPlez et al.l 1 19931 : lAbia et ail 
1999), this line corresponding to the transition 2p-3s. It 
has also been suggested that this line might be visible 
in other Li-rich objects such as chemically peculiar stars 
(jPolosukhina et aT1ll999ft . though we are not aware of any 
such detections. 

Previous non- LTE calculation s of Li line formation 
in co ol stars by ICarlsson et al.l 



( 1994ft and iLind et al.l 
(|2009ft used the rate coefficients of iPar Id (j 1971ft for the de- 
scription of the electron-im pact e x citatio n in their model 
atom and the recipe of [ Allen! (|1976ft for ionization. 
ISteenbock fc Holwegerl (|1984 ) also used the Park data, but 
replaced the value for the 2s-2p transition with the van 
Regemorter value. A similar semi-empirical recipe for the 
ionization is also used. Thus, the three works that are of 
most interest for comparison use essentially the same elec- 
tron colli sion data. In th is work we used the same model 
atom as lLind et al.l (|2009ft . changing only the electronic col- 
lision data; that is, the rate coefficients of Tables [3] and [5] 
replace the old data for excitation an d ioni z ation by elec- 
tron impact. The collision rates from iParkl ()1971l ) for the 
transitions involving high-lying levels not covered by our 
calculations were retained. We calculated non-LTE correc- 
tions for a range of stellar atmospheric models using the 
old and new model atoms. Some example comparisons of 
the results are shown in Fig. @] We found the differences 
between the two model atoms to be < 0.01 dex for solar- 
type stars. In warm stars, T e g > 7000 K, the differences are 
slightly larger but still < 0.02 dex. 

The behaviour of the abundances measured using the 
6104 A and 8126 A spectral lines was also studied. The 
transition 2p-3s, corresponding to the 8126 A line, has the 
largest differen ce in the ra te coefficients when compared 
with those from lParkl (|l97ll ). For all the atmospheric model 
cases shown in Fig. [H the Li abundance derived from this 
line shows very similar sensitivity to the new electron colli- 
sion data. The effects of the new electron collision data on 
the abundance measured from the 6104 A line are consider- 
ably less in the solar case: only about 10% of the difference 
shown by the other two lines. In the other cases the influ- 



ence of the electron collision data on the abundance mea- 
sured from the 6104 A spectral line is again comparable to 
the influence on the 6708 and 8126 A lines. 



4. Conclusions 

A common problem for calculating non-LTE abundance 
corrections is that they can be subject to errors from un- 
certainties in the underlying atomic data. We have shown 
that electron collision data from a number of different cal- 
culations based on quite different methods are in good or 
even excellent agreement. Data calculated with advanced 
close-coupling techniques provides data in excellent agree- 
ment, giving rate coefficients for excitation processes within 
a factor of 2 for the temperatures of most interest in cool 
stellar atmospheres, T ~ 5000 - 8000 K. This, and the good 
agreement with experi mental results where available (see 
Schweinz er et al.lll999ft . suggests the data have an uncer- 
tainty of similar magnitude. The semi-empirical data for 
both excitation and ionization that have been commonly 
used in non-LTE calculations for Li up to now, agree rea- 
sonably with the modern data, generally agreeing within 
a factor of 6, though with a few significant outliers. This 
suggests that the estimates of the uncertainties of a factor 
of 2 for these data are somewhat too low, though the data 
for the important 2s-2p transition does indeed agree within 
a factor a two. 

We find that these uncertainties are not very impor- 
tant when applied to Li I in cool stars, and the differences 
typically result in uncertainties of less than 0.01 dex. The 
strongest effects are found in F dwarfs and extremely Li-rich 
stars. These uncertainties are certainly negligible compared 
to those arising fr om other sour ces, especially the atmo- 
spheric modelling (Asplund 2005). Thus, we conclude that 
the existing electron collision data, or the data presented in 
this paper, when coupled with data for the charge exchange 
process Li(3s) + H ^ Li+ + IT now available, are accurate 
enough to reliably model non-LTE formation of Li lines in 
stellar atmospheres at the 0.01 dex (~ 2%) level. 
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Table 5. Rate coefficients for collisional ionization between electrons and Li I atoms for different temperatures. 



Bound 






Temperature (K) 






State 


1 000 


2 000 


5 000 8 000 12 000 


20 000 


50 000 



2s( 2 S c ) 
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2.73(- 


-09) 
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Notes. The rate coefficients (<rv), in cm 3 s 1 , are calculated from the ionization cross sections from ILennon et al.l (|1988l ) and 
ISchweinzer et al.l (|1999T l (see text). The values should be read as: a(b) := a x 10 b 
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Fig. 4. Non-LTE Li abundance c orrections fo r the 6707 A resonance line using the model atom as lLind et al.l (|2009[ ) and 
replacing the e+Li collisions from lParkl (|l97lh (dashed line) with the new calculations (solid line). The dotted line shows 
the zero value. In the lower left corner of each plot are the stellar parameters of the model used: T e ff/log g/[Fe/H}. 
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